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Abstract 

Within  the  framework  of  density  functional  theory,  we  investigate  the  structural  properties,  electronic  structure,  and  the  chemical  stability  of 
magnesium  alanate  (Mg(AlH4)2),  a  complex  hydride  for  hydrogen  storage,  using  the  projector  augmented  plane  wave  method.  The  calculated 
structural  properties  of  magnesium  alanate  compare  well  with  the  recent  experiment.  A  detailed  study  of  the  electronic  structure  and  charge 
density  redistribution  reveals  the  orbital  hybridization  and  characteristics  of  bonding  orbitals  within  the  complex  hydride.  We  found  an  ionic 
covalent  bonding  between  A1  and  H  in  the  (A1H4)~  anionic  complexes  embedded  in  the  matrix  of  Mg2+  cations.  The  calculated  reaction  energy  of 
Mg(AlH4)2  is  in  good  agreement  with  the  experimentally  determined  value. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  use  of  hydrogen-based  energy  in  practical  applications 
such  as  fuel  cell  vehicles  requires  the  development  of  safe  and 
efficient  hydrogen  storage  technology.  At  present,  there  are  three 
generic  routes  for  the  solid  state  storage  of  hydrogen  [1,2]; 
(i)  physisorption  as  in  many  porous  carbon  and  zeolite  mate¬ 
rials  with  high  specific  surface  area;  (ii)  hydrogen  interaction  in 
metal  hydrides;  (iii)  chemical  reaction  such  as  in  complex  metal 
hydrides.  Group  1,  2,  and  3  light  metals  (e.g.  Li,  Mg,  B,  and  Al) 
give  rise  to  a  large  variety  of  metal-hydrogen  complexes,  which 
are  especially  interesting  and  have  been  attracted  considerable 
attention  because  of  their  potential  use  as  lightweight  materi¬ 
als  for  reversible  hydrogen  storage  [1  ].  Compared  with  LiAlH4 
and  NaAlH4  (with  10.6  and  7.5  wt.%  theoretical  hydrogen  con¬ 
tent,  respectively),  magnesium  alanate  (Mg(AlH4)2),  is  another 
potential  hydrogen  storage  candidate  with  a  hydrogen  content  of 
9.3  wt.%  [3].  This  compound  decomposes  in  a  two-step  reaction 
[3]: 

Mg(AlH4)2  -*  MgH2  +  2A1  +  3H2  (1) 

MgH2  —>■  xa  +  y(3  +  H2  (2) 

where  a  and  |3  in  reaction  (2)  are  Al-Mg  solid  solutions  (such 
as  ALMgo).  Reaction  (1)  has  been  found  to  occur  at  temper- 
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atures  of  135-163  °C  [4-7]  and  reproduce  7.0  wt.%  released 
hydrogen.  The  second  reaction  step  has  been  reported  that 
it  companies  with  2.3  wt.%  released  hydrogen  and  occurs  at 
270-3 10 °C  [5,8],  However,  the  formation  of  stable  MgH2  as 
an  intermediate  or  AI2Mg2  as  an  end  product  might  negatively 
influence  the  reabsorption  of  hydrogen  when  starting  from  the 
respective  dehydrogenated  state  of  the  material  [7].  In  order  to 
improve  the  performance  of  Mg(AlH4)2  involved  in  reversible 
hydrogen  absorption/desorption  processes,  understanding  the 
fundamental  physical  properties  of  this  complex  hydride  is  of 
key  importance  and  great  interest. 

Although  theoretical  investigations  using  the  first-principles 
calculations  are  so  common  for  alkali  complex  hydrides 
[2,9-1 1],  only  few  studies  on  Mg(AlH4)2  has  been  reported,  in 
which  one  [12]  a  variety  of  Mg(AlH4)2  clusters  were  served 
as  models  of  the  bulk  structure  and  it  is  mainly  concerned 
with  the  structural  properties.  The  reaction  enthalpy  of  hydro¬ 
gen  desorption  for  Mg(AlH4)2  is  estimated  as  41  kJ  moL  1  per 
H2  molecule  based  on  the  cluster  model  [12]  or  as  46  kJ  mol-1 
per  H2  molecule  from  ab  initio  vibration  properties  calculations 
based  on  density-functional  perturbation  theory  [13],  respec¬ 
tively,  which  are  both  much  larger  than  the  value  of  0  kJ  mol- 1 
suggested  in  the  review  study  [14]  in  light  of  the  experimental 
measures  [14,15]. 

Therefore,  we  employ  first-principles  approaches  based  on 
density  functional  theory  and  the  projector  augmented  plane 
wave  method  to  conduct  a  comprehensive  study  of  the  structural, 
electronic,  and  stability  properties  of  Mg(AlH4)2.  The  system- 
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atic  analysis  of  the  electronic  properties  of  Mg(AlH4)2  suggests 
a  covalent  type  of  bonding  between  A1  and  H  within  the  complex, 
but  with  a  large  ionicity.  This  is  the  generic  bonding  characteris¬ 
tics  in  alkali  complex  hydrides  [9] .  From  the  static  total  energy 
calculations,  i.e.  the  zero-point  energies  are  not  included,  we 
have  obtained  that  a  reasonable  reaction  enthalpy  of  hydrogen 
desorption  for  Mg(AlH4)2  is  — 2.22kJ  mol-1  which  is  in  good 
agreement  with  the  experimental  value. 

2.  Computational  methodology 

First-principles  calculations  on  Mg(AlH4)2  were  performed 
using  the  Vienna  ab  initio  simulation  package  (VASP)  [16,17], 
which  is  based  on  the  density  functional  theory  (DFT)  [18].  The 
generalized  gradient  approximation  (GGA)  of  Perdew  et  al.  [19] 
known  as  GGA-PW9 1  was  used  to  treat  the  electronic  exchange- 
correlation  energy.  Electron-ion  interaction  was  represented  by 
the  projector  augmented  wave  (PAW)  method  [20]  with  plane 
waves  up  to  an  energy  of  550  eV.  The  reference  configurations 
for  valence  electrons  were  3s23p°  for  Mg,  3s“3p 1  for  Al,  and 
Is1  for  H,  respectively. 

To  examine  the  transferability  of  the  pseudopotentials,  the 
calculations  have  been  performed  on  the  related  solid  (Mg,  Al, 
and  MgH2)  and  a  H2  molecule.  We  choose  the  rutile-type  a- 
MgFE  for  simplicity.  The  same  energy  cutoff  with  GGA-PW91 
is  used.  For  solids,  the  k-point  grids  for  Brillouin  zone  inte¬ 
gration  are  generated  according  to  the  Monkhorst-Pack  scheme 
[21]  so  as  to  make  the  edge  lengths  of  the  grid  elements  closer 
to  the  target  value  of  0.08  bohr-1  (1  bohr  =  0.529  A)  as  possi¬ 
ble.  The  same  conditions  are  applied  for  the  calculations  on 
Mg(AlH4)2.  The  resulting  k-point  grid  is  9  x  9  x  9  (349  irre¬ 
ducible  points)  for  hexagonal  MgfAlFUh.  We  have  carefully 
checked  that  these  conditions  give  good  convergence  of  the  total 
energy  within  1  me V  atom- 1 .  The  calculation  for  an  isolated 
H2  molecule  is  carried  out  using  a  cubic  supercell  with  size  of 
15  A  x  15  Ax  15A  and  with  the  single  /’  point  for  the  k-point 
sampling.  The  results  are  given  in  Table  1.  The  lattice  constants 
(interatomic  distance),  atomic  positions  and  cohesive  energies 

Table  1 

Results  for  hcp-Mg,  fcc-Al,  H2  molecule  and  rutile-type  a-MgH2,  including  the 
lattice  constants  a  and  c  (in  A),  cohesive  energy  Ecoh  (in  eV  per  formula  unit), 
interatomic  distance  d  (in  A),  and  positions  of  H  (wh)  atoms 


Material 

Property 

This  work 

Experiment 

Reference 

hcp-Mg 

a 

3.18 

3.21 

[27] 

c 

5.24 

5.21 

Ecoh 

1.48 

1.51 

fcc-Al 

a 

4.04 

4.05 

[27] 

Ecoh 

3.49 

3.39 

h2 

cl 

0.75 

0.741 

[28] 

Ecoh 

4.56 

4.74 

a-MgH2 

a 

4.495 

4.501 

[29] 

c 

3.006 

3.010 

uu 

0.304 

0.304 

Ecoh 

6.70 

- 

Note  that  the  zero-point  energies  are  not  included  in  Ecoh- 


Fig.  1.  Crystal  structure  of  hexagonal  Mg(AlH4)2. 


of  the  related  solid  and  molecule  are  in  good  agreement  with  the 
available  experimental  data. 

3.  Results  and  discussion 

3.1.  Crystal  structure 

Mg(AlH4)2  crystallizes  in  a  hexagonal  structure  having  the 
P3m  1  space  group.  The  structure  was  confirmed  by  combined 
powder  synchrotron  X-ray  and  neutron  diffraction  [3,7].  The 
ball-and-stick  model  for  crystal  structure  of  Mg(AlH4)2  is  illus¬ 
trated  in  Fig.  1 .  The  structure  consists  of  isolated  and  slightly  dis¬ 
torted  AIH4-  tetrahedra  that  are  connected  via  six-coordinated 
Mg  atoms  in  a  distorted  octahedral  geometry,  resulting  in  a 
sheet-like  structure  along  the  crystallographic  c-axis  [3].  Start¬ 
ing  from  experimentally  available  parameters,  the  lattice  vectors 
of  primitive  cell  and  atomic  coordinates  were  relaxed  at  a  series 
of  fixed  volumes  by  optimizing  both  the  forces  and  stresses.  The 
residual  force  and  stress  in  the  equilibrium  geometry  are  of  the 
order  of  1  meV  A-1  and  10-3  GPa,  respectively.  The  obtained 
energies  were  fitted  with  a  Mumaghan  equation  of  state  [22]  to 
give  the  equilibrium  volume  and  the  minimum  energy.  In  order  to 
further  asses  the  effects  of  the  exchange-correlation  approxima¬ 
tion  on  the  structural  properties  the  local  density  approximation 
(LDA)  [23]  and  GGA  of  Perdew-Burke-Ernzerhof  (PBE) 
parameterization  [24]  were  also  considered  in  addition 
to  GGA-PW91.  The  final  calculated  cell  parameters  and 
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Table  2 

Comparison  between  calculated  and  experimental  lattice  constants,  internal 
parameters,  and  bond  lengths  for  hexagonal  Mg(AlH4)2  (space  group  P3ml) 
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This  work 

LDA 

GGA-PW91  GGA-PBE 

Experimental  [3]  (8  K) 

Lattice  constants  (A) 

a 

5.133 

5.232 

5.245 

5.208 

c 

5.594 

6.025 

6.044 

5.839 

Volume  (A3) 

127.64 

142.86 

144.01 

137.16 

Internal  parameters 

Af 

0.688 

0.706 

0.706 

0.6991 

HI, 

0.340 

0.440 

0.441 

0.4242 

H2, 

0.165 

0.168 

0.168 

0.1671 

H2C 

0.796 

0.812 

0.812 

0.8105 

Bond  lengths  (A) 

Mg-H2 

1.86 

1.90 

1.91 

1.870 

Mg-Al 

3.44 

3.50 

3.51 

3.482 

Al-Hl 

1.61 

1.60 

1.60 

1.606 

A1-H2 

1.62 

1.63 

1.63 

1.634 

The  Wyckoff  positions  of  Mg,  Al,  HI,  and  H2  atoms  are  1  a  (0,  0,  0),  2d  (1/3, 
2/3,  z ),  2d  (1/3,  2/3,  z ),  and  6 i  (x,  —  x,  z ),  respectively,  x  and  z  are  the  fractional 
coordinates  of  the  nonequivalent  sites  of  Al  and  H  atoms  in  Mg(AlH4)2- 


Fig.  2.  Energy  band  structure  of  Mg(AlH4)2.  The  energy  zero  is  set  at  the  top 
of  the  valence  band  and  shown  as  dotted  line. 


bond  lengths  are  given  in  Table  2,  along  with  experimental 
values. 

As  seen  in  Table  2,  GGA-PW91  and  GGA-PBE  consistently 
yields  the  larger  volumes  (4%  and  5%,  respectively)  than  exper¬ 
iment,  which  are  typical  for  the  GGA  approximation  to  DFT. 
While  LDA  gives  a  smaller  (7%)  volume  than  that  of  experiment. 
We  notice  that  GGA-PW9 1  also  yields  a  better  cell  parameters  of 
Mg(AlH4)2,  including  the  lattice  constants,  positions  of  H  atoms 
and  the  bond  lengths.  For  example,  the  reported  Mg(AlH4)2  lat¬ 
tice  constants  (at  8K)  [3]  are  5.208  and  5.839  A,  respectively, 
while  we  obtain  5.232  and  6.025  A.  Our  calculated  Al  position 
(1/3,  2/3,  0.706)  and  H  positions  (1/3,  2/3,  0.440)  and  (0.168, 
—0.168,  0.812)  are  also  in  good  agreement  with  the  reported 
values  [3]. 


aluminum  exists  in  the  anion  complex  AIH4- .  This  characteris¬ 
tic  bonding  of  complex  AIH4  has  also  been  reported  in  NaAlFL; 
compound  [9]. 

In  order  to  further  analyze  the  bonding  in  Mg(AlH4)2,  the 
charge  distribution  and  charge  transfer  in  the  AIH4  complexes 
are  also  examined.  The  total  charge  density  p(r)  and  the  dif¬ 
ference  charge  densities  Ap(r)  in  the  plane  containing  Mg,  Al, 
and  H  atoms,  are  shown  in  Fig.  4.  The  difference  charge  density 
Ap(r)  [9]  is  defined  as  the  difference  between  the  total  charge 
density  of  the  solid  and  a  superposition  of  atomic  charge  densi¬ 
ties  with  the  same  spatial  coordinates  as  in  the  solid.  In  Fig.  4,  the 
total  charge  density  plot  shows  that  the  buildup  of  the  electron 
density  is  mainly  around  the  positions  of  H  atoms  and  exhibits 
a  slightly  directional  distribution  toward  Al  atoms,  while  the 
difference  density  plot  shows  the  charge  transfer  from  Al  to  H 


3.2.  Electronic  structure 

The  energy  band  structure  and  electronic  density  of  states 
(DOS)  ofMg(AlH4)2  are  calculated  and  shown  in  Figs.  2  and  3, 
respectively.  Both  of  these  figures  exhibit  a  rather  large  energy 
gap  of  about  4.3  eV  between  valence  and  conduction  bands,  indi¬ 
cating  that  Mg(AlH4)2  is  a  wide-gap  insulator.  The  eight  valence 
bands  in  Fig.  2  have  a  width  of  about  6  eV  and  are  split  into  two 
disjointed  groups  consisting  of  two  and  six  bands,  respectively. 
These  two  groups  do  not  overlap  in  energy.  From  these  figures  it 
can  be  seen  that  the  first  two  valence  bands  from  —6  to  —4  eV  are 
derived  from  FI  s  and  Al  s  orbitals,  while  the  next  six  bands  in  the 
higher-energy  group  between  —3  and  0  eV  are  mainly  from  H  s 
and  Al  p  states.  This  suggests  that  there  is  a  strong  hybridization 
between  H  and  Al  orbitals.  The  conduction  bands  mainly  by  Mg 
(and  to  lesser  extent  Al)  empty  states.  Mg  has  hardly  any  pro¬ 
jection  in  the  occupied  states,  while  H  has  hardly  any  projection 
in  the  range  of  the  empty  states  shown  in  Fig.  3.  These  features 
indicate  that  a  highly  ionic  covalent  bonding  of  hydrogen  and 
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Fig.  3.  Electronic  density  of  states  (DOS)  of  Mg(AlH4)2,  including  the  total 
and  projections  DOS.  The  energy  zero  is  set  at  the  top  of  the  valence  band  and 
shown  as  dotted  line. 
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Fig.  4.  Total  and  difference  charge  density  plots  of  Mg(AlH4)2  in  the  (110) 
plane,  (a)  The  total  charge  density  represented  by  contours  in  an  interval  of 
0.09  electrons  A~3  from  0  to  1 .87  electrons  A  \  (b)  Charge  deficiency  in  the 
difference  plot  is  plotted  in  dashed  lines,  while  the  density  increase  near  the 
hydrogen  atoms  is  plotted  in  solid  lines  with  an  increment  of  0.06  electrons  A-3 
from  —0.108  to  0.63  electrons  A^3 . 

leaving  aluminum  positively  and  hydrogen  negatively  charged. 
These  features  support  a  covalent  type  of  bonding  but  with  a 
large  ionicity  in  AIH4  complexes. 

3.3.  Formation  energy 

The  heat  of  formation  is  the  most  fundamental  and  important 
quantity  for  the  hydrogen  storage  materials.  For  MgiAlFL^,  the 
cohesive  energy  A’coh  is  calculated  as  27.36  eV  per  formula  unit. 
Using  this  value  and  the  energies  given  in  Table  1,  the  heat  of 
formation  for  Mg(AlH4)2  is  calculated  according  to  the  reaction 
as 

Mg  +  2A1  +  4H2  -+  Mg(AlH4)2  (3) 

If  only  the  total  energies  are  considered,  the  heat  of  forma¬ 
tion  for  Mg(AlH4)2  is  estimated  as  +0.63  eV  per  formula  unit 
(i.e.  — 61.07kJ mol-1).  In  comparison,  the  heat  of  formation 
for  a-MgH2  is  obtained  — 63.36kJ  mol-1,  which  is  in  excellent 
agreement  with  the  previous  theoretical  study  (— 64kJmol-1 
without  zero-point  energy  contributions)  [25]  and  in  good  agree¬ 
ment  with  the  experiment  value  (—75.2,  —67.0,  —74.4,  and 
—65.8  kJ  mol-1)  [26].  Following  Eq.  (1)  which  is  more  closely 
related  to  hydriding  and  dehydriding  reactions,  it  is  possible 


to  determine  the  reaction  enthalpy  of  hydrogen  desorption  for 
Mg(AlH4)2.  Therefore  the  reaction  enthalpy  of  hydrogen  des¬ 
orption  for  Mg(AlH4)2  is  predicted  as  — 2.22kJ  mol-1  without 
the  zero-point  energies,  which  is  also  very  close  to  the  value 
of  Oklmol-1  H2  reported  in  experiment  [14,15].  It  indicates 
that  the  thermodynamic  properties  of  Mg(AlH4)2  appear  the  be 
unfavorable  for  reversible  hydrogen  storage  under  practical  con¬ 
ditions  [14]  due  to  stable  MgH2  as  an  intermediate. 

4.  Conclusions 

The  structural,  electronic  and  stability  properties  of 
Mg(AlH4)2  have  been  calculated  and  analyzed.  Our  calculated 
structural  parameters  including  the  lattice  constants  and  internal 
parameters  are  in  good  agreement  with  the  experiment  results. 
Especially,  the  hydrogen  positions  in  Mg(AlH4)2  are  very  close 
to  the  deuterium  positions  obtained  in  the  neutron  diffraction 
measurement  at  8K.  The  volume  of  the  unit  cell  calculated 
with  the  LDA  approximations  is  about  7%  smaller  than  the 
experimental  value.  The  best  calculated  volume  for  Mg(AlH4)2 
is  obtained  with  the  GGA-PW91  approach.  Detailed  calcula¬ 
tions  of  electronic  structure  suggest  that  the  complex  hydride 
Mg(AlH4)2  is  an  insulator  with  GGA  energy  band  gap  of  4.3  eV. 
The  total  and  projected  electronic  DOS  combined  with  the 
charge  density  plots  indicate  significant  charge  transfer  lead¬ 
ing  to  Mg2+  ions  and  negatively  charged  AIH4  complex.  The 
bonding  within  the  AIH4  complex  involves  the  hybridization 
between  A1  and  H  orbitals,  but  shows  a  strong  ionic  charac¬ 
ter.  Total  energy  calculations  are  used  to  estimate  the  enthalpy 
change  in  the  chemical  reactions  associated  with  the  hydra¬ 
tion/dehydration  processes.  The  calculated  values  are  very  close 
to  the  measured  values  in  experiments. 
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